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ABSTRACT 



Context. Lyo' is a key diagnostic for numerous observations of distant star-forming galaxies. It's interpretation requires, however, 
detailed radiation transfer models. 

Aims. We provide an extensive grid of 3D radiation transfer models simulating the Lyce and UV continuum radiation transfer in the 
interstellar medium of star-forming galaxies. 

Methods. We have improved our Monte Carlo McLya code, and have used it to compute a grid of 6240 radiation transfer models for 
homogeneous spherical shells containing Hi and dust surrounding a central source. The simulations cover a wide range of parameter 
space. We present the detailed predictions from our models including in particular the Lya escape fraction f^^^, the continuum atten- 
uation, and detailed Lyo- line profiles. 

Results. The Lyo' escape fraction is shown to depend strongly on dust content, but also on other parameters (Hi column density and 
radial velocity). The predicted line profiles show a great diversity of morphologies ranging from broad absorption lines to emission 
lines with complex features. The results from our simulations are distributed in electronic format. 

Conclusions. Our models should be of use for the interpretation of observations from distant galaxies, for other simulations, and 
should also serve as an important base for comparison for future, more refined, radiation transfer models. 

Key words. Galaxies: starburst - Galaxies: ISM ~ Galaxies: high-redshift - Ultraviolet: galaxies ~ Radiative transfer 
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1. Introduction 

Stimulated by numerous observations of galaxies and the 
intergalactic medium (IGM) in the high redshift Universe, 
and made possible by significant growth of computational 
power, several groups have recently d eveloped numerical codes 
treatin g the transf er of Lya radiation (Zheng & Miralda- Escude 
20021: lAhn etall \200^ iDiik stra et al. 2006; Verhamme et al 
2006 ; Ij asitsiomi: 2006 ; Laursen & Sommer; LarsenI 120071 
Barnes & Haehnefl 2009 . e.g.). Most of the codes use the 
Monte Carlo (MC) method to sample the radiation transfer 
The main advantages of this technique is its ease in imple- 
mentation, treatment of arbitrary geometries, and inclusion 
of different physical scattering and absorption processes. The 
codes have been applied to a variety of astrophysical problems, 
including e.g. cosmic reionisation. Lyes' fluorescence in the 
IGM, predictions of Lya emission from distant galaxies, and 
the interpretation of observed Lya line profil es from starburst 
galaxies (Cantaluoo et al. 2005, 2008; Sem eUn et all l2007t 
Baeketal. 2009;' Pierleoni et al. 2009; Barnes & Haehneltl 
2010; KollmeieLet al. 2010; Laursen et aL ,201 k Zheng et alJ 
20lil). 



Our group has recently developed a general 3D radi- 
ation transfer code including Lya and UV continuum ra- 
diation including scattering and absorpt ion by neutral hy- 
drogen and dust dVerhamme et al.l l2006l) . So far, our code 



has mostly been applied to understand the diversity of ob- 
served Lyo' line profiles and to quantitatively exploit the Lya 
profile in Lyo' emitters (LAE) and Lyman Break Galaxies 
(LBG) (seel^h aerer & Verham mel2008l;]Ve i-hamme et al. ll2008l 
iDessauges-Zavadskv et al.i i2010l : Vanzella et al. 2010) and for 
the modeling of local starburst galaxies (cf. Atek et al. 20091^. 
In particular, we have shown that the observed diversity of line 
profiles reaching from broad absorption in some LBGs to asym- 
metric emission in LAEs, can be understood by variations of 
the Hi column density and dust content in spherically expanding 
shells. Furthermore our models also place interesting constraints 
on the intrinsic Lya emission properties, hence on age and star- 
formation history. For example, for LBGs we found no need for 
particular ages, and we showed that their intrinsic Lyo- equiva- 
lent widths are compatible with constant star-formation over sev- 
eral hundred Myr, as also indicated by their broad-band SEDs. 
Last, but not least, our simulations have clarified the relation be- 
tween ob served velocity shif ts of Lya and ISM expansion velo c- 
ities (see ISchaerer & Verham me 2008; f^rhamme et al.ll2008l ). 

These studies have been based on a relatively small num- 
ber of simulations and tailored models to fit observations of less 
than twenty galaxies, albeit of different types. To explore a larger 
range of the parameter space, and to do this in a systematic way, 
allowing us also to examine possible degeneracies in Lya line 
profile fits, we have computed a large grid of radiation trans- 
fer models for homogeneous, spherically expanding shells. Our 
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simulations cover a wide range of the 4-dimensional parameter 
space given by the expansion velocity, Hi column density, ve- 
locity dispersion (b), and dust optical depth. From this, spectra 
can be predicted for any arbitrary input spectrum, containing e.g. 
a UV continuum plus intrinsic Lya emission from the central 
source. 

To make these simulations feasible, we have parallelised our 
Monte Carlo code, named McLya. Furthermore we have added 
some new physics in the code, such as scattering by Deuterium 
and the recoil effect, and we have made other improvements (in 
particular on scattering by dust). Our standard code works on 
cartesian grids; a version using nested-grids, provided by AMR 
(Adaptive Mesh Refinement) codes has recently been developed 
by Verhamme et al. (in preparation). In addition we have devel- 
oped an automated tool to fit observed Lya line profiles, making 
use of our grid of McLya simulations. 

The description of the grid of UV and Lya radiation transfer 
simulations and some results derived from it, are the main ob- 
jective of this paper. At the same time the model grid described 
here will be made available electronically. 

Predictions from our mo del grid have alread y been used in 
several papers. For example, lAtek et all (l2009al) compared pre- 
dicted LycK escape fraction from our grid t o measurem e nts of 
this q uantity in z ~ 0.2-0.4 galaxies, and iHaves et"^ (1201 Ol 
1201 li) have compared the same quantity to measure ment s atred- 
shift z ~ 2 and hig her In .Dessauges-Zavadskv et al.l (l2010h : 
IVanzellaetan(l2O10() . we have used our fitting tool to analyse 
two individual LBGs at redshift z = 2.8 and z = 5.8 respec- 
tively. The predictions from our model grid have recently been 
included in the sem i-analytical galaxy models of the Lyon group 
dGarel et al.ll201 lb . Our predictions are also being used to pre- 
pare observations with future instruments, such as MUSE for 
the VLT (Garel et al., in preparation). 

Our paper is structured as follows. In Sect. |2]we describe 
the main improvements implemented in our code. In Sect.|3]we 
summarise the input parameters and other issues for the compu- 
tation of our models grid. Predictions for the Lya escape fraction 
and the strength of Lya absorption are shown in Sect. |4] In Sect. 
|5] we illustrate some line profiles predicted from our model grid. 
In Sect. |6] we summarise our main conclusions. 

2. Update of the MCLya code 

We have developed an improve d version of the M onte Carlo ra- 
diation transfer code McLya of i Verhamme et al.l (|2006) includ- 
ing the detailed physics of Lya line and UV continuum transfer, 
dust scattering, and dust absorption for arbitrary 3D geometries 
and velocity fields. 

The following improvements have been included: 

- Angular redistribution functions taking quantum mechani- 
cal results for Lya into account. More precisely we dis- 
tinguish two different phase f unctions for scattering i n the 
core and the wing, following iDiikstra & Loebl (l2008l) . For 
the distinction between these two regimes we also use their 
value of Xcrit = 0.2, where x is the frequency shift from line 
center expressed in Doppler units. In the wing, the phase 
function is identical to the one of dipolar (Rayleigh) scat- 
te ring, adopted previously in our code for the computations 
of I Verhamme et all (1200 8). Although the differences in the 
phase functions are important for a proper treatment of po- 
larised radiation (cf. IDiikstra & Loebl2008l) . we have not no- 
ticed changes in our predictions. 



- Freq uency changes of Lyg photons due to the recoil effect 
(e.g. lZheng & Miralda-Escudel200 2'). Before, scattering was 
considered coherent in the atom's frame. As well known, the 
recoil effect can lead to differences for very low t emper - 
atures, as shown e.g. by IZheng & Miralda-Escudil (|2002|) : 
iTasitsiornil (120061) . 

- The presence of Deuterium. Following the suggestion of 
Ipijkstra et al. (2006), Deuterium has been included, assum- 
ing a canonical abundance of DjH = 3 x 10"^. It's effect 
may be visible for cases with a static ISM in the blue part 
of the line (cf. Fig. 3 of IDiikstra et al.ll2006h . For expand- 
ing geometries, such as the ones discussed below, it's effect 
becomes, however, very small or invisible. 

- Anisotropic dust scattering using the Henyey-Greenstein 
phase function. We adopt the following values for the dust 
albedo Qa - 0.46 and the scattering phase function asym- 
metry g - (cosO) = 0.77 at Lya following I Witt & GordonI 
(2000). In our previous computations we adopted Qa = 0.5 
and g = (isotropic dust scattering). Among the changes 
made, this improvement has the largest impact on the results 
since it leads to a more pronounced fore-aft symmetry for 
dust scattering, which in turn affects the escape fraction of 
photons both in the UV continuum and in the Lya line. For 
the spherical shell geometry considered below, this leads in 
general to somewhat higher escape fractions, i.e. less effi- 
cient dust attenuation. 

- Finally, the code has been parallelised using OpenMPI, for 
efficient use on supercomputers. Given the use of the Monte 
Carlo method, parallelisation is basically trivial, and near- 
perfect scaling properties are achieved. Typically, we have 
run our simulations on 100-300 CPU cores in parallel. 

Th e remaining input physics is described in lVerhamme et all 
( |2006|) . We do not rely on approximations to accelerate the 
Lyo' transfer Although we have experimented with some meth- 
ods (e.g. those in Verhamme et al. 2006; Laursen et al. 2009al), 
we have noted some differences in the resulting line profiles. 
Therefore, to avoid possible inaccuracies in our calculations ex- 
ploring a large parameter space, we prefer to use a correct, "brute 
force" approach. 

3. Grid of MCLya models 

For simplicity, and to provide a basis for other future studies, 
all simulations carried out subsequently assume a homogeneous 
and co-spatial distribution of neutral hydrogen and dust with a 
constant density and temperature. The effects of clumpy media 
are examined in an upcoming publication (Duval et al., in prepa- 
ration). The input parameters of the code are the Hi geometry 
and velocity field, the microscopic Hi velocity distribution, the 
spatial location and distribution of the UV continuum and line 
emission source(s), and the dust-to-gas ratio. 

For the current model grid we consider the following ge- 
ometry: spherically symmetric shells with a central source UV 
source. This case is described by 4 parameters: (i) the radial ex- 
pansion velocity of the shell, Vexp, (H) the Hi column density 
towards the source, A^hi, (Hi) the microscopic Hi velocity distri- 
bution described by the Doppler parameter b, and (iv) the dust 
absorption optical d epth Tg, which ex presses the dust-to-gas ra- 
tio. As discussed by I Verhamme et all (2006) t,, is related to the 
usual color excess E{B -V)hy E{B - V) ~ (0.06.. . 0.1 1)t„. For 
the Calzetti et al. (2000) law one has E{B-V) « 0.1 r„. In short, 
the present model grid is described by 4 input parameters Vexp, 
A^Hi, b, and r^. By symmetry of the radiation transfer problem 
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Table 1. Input parameter values adopted for the grid of radiation transfer models of spherically expanding shells. All combinations 

^exp ' 



of parameters have been computed, corresponding to 6240 models in total. The units are the following: Vgxp and b in km s \ A^hi in 



cm ^. is dimensionless. By symmetry of the problem, the models also describe spherical infall, in which case the infall velocity 

is V = -Vexp. 

Parameter Values n 



Vexp 0.,20., 50., 100., 150., 200., 250., 300., 400., 500., 600., 700. 12 

log(iVHi) 16., 18., 18.5, 19.,19.3,19.6,19.9,20.2,20.5,20.8,21. 1,21.4,21.7 13 

b 10., 20., 40., 80., 160. 5 

Ta 0., 0.2,0.5, 1., 1.5,2., 3., 4. 8 



(cf. iNeufeldl [19901; I Verhamme etaU 120061 iDiikstra et all 120061) 

our grid calculations also apply to the case of spherical infall. 
In other words models with a radial expansion velocity Vexp also 
describe the case of spherical infall with v - -Vexp. 

Thanks to the parallelisation of the code, the computation of 
a large grid, covering a wide 4-dimensional parameter space, has 
been possible. The values adopted for the input parameters are 
listed in TablelT] Simulations have been run for all combinations, 
yielding in total 6240 models. 

For each parameter set, a full Monte Carlo simulation is run 
with 1000 photons per input frequency bin. The radiation trans- 
fer calculations cover a sufficiently broad spectral range (here 
typically from -6000 to +6000 km s"', in bins of 20 km s"' 
for the input) to reach the continuum for most simulations (ex- 
cept for s ome of the highes t column density simulations). As de- 
scribed in lVerhamme et al](f2006 ). our MC simulations are com- 
puted for a flat input spectrum, keeping track of the necessary 
information to recompute a posteriori simulations for arbitrary 
input spectra. The total computing time required for the entire 
grid amounted to ^ 25 CPU year. 

In practice, the following quantities are stored for each emer- 
gent photon: input and output frequency, position and direc- 
tion of escape. From this we can in particular compute the pre- 
dicted emergent spectrum (integrated or spatially resolved spec- 
tra, spectral maps) and the escape fraction for all frequency bins. 
The detailed model results (approx. 60 GB) allowing one to com- 
pute emergent spectra are available on request from the first au- 
thofl Derived quantities (the Lya escape fraction computed for 
different FWHM, the continuum escape fraction, and the equiv- 
alent width of the Lya absorption assuming a flat continuum, cf. 
below) for all models are provided in Table |2] 

4. Predictions for tlie Lya escape fraction and 
related quantities 

The total Lya escape fraction, /esc, is defined as the ratio be- 
tween the number of Lya line photons emitted, A^, by the UV 
source and the number of these photons, A'esc, emerging from the 
simulation box. The monochromatic escape fraction is given by 
/esc(/l) = Nesc{^)/N(A), where A is the "input" wavelength, i.e. 
that of the photons emitted from the source before being altered 
by radiation transfer The total Lya escape fraction is 

^'"'^ ~ jN(A)dA ' 

where the integration is carried out over the entire line profile 
of the "input" spectral line. In principle /esc therefore depends 
on the emission line profile of the source. In practice, however, 
/esc is close to the monochromatic escape fraction at line center 

' See,http : //obswww . unige . ch/sf r tOr contact DS. 
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Fig. 1. Predicted Lya escape fractions from our 6240 models 
as a function of Eb-v- Yellow (red) symbols show the simula- 
tions with dust-to-gas ratios within a factor 10 (2) of the Galactic 
value. Among the models with dust-to-gas ratios within a factor 
2 of the Galactic value, those with Vexp=200 km s"' are shown 
by blue symbols. Black symbols show all the remaining models. 
The blue dashe d line shows an att enuation with = 12. corre- 
sponding to the 'Calzett i et al.l ( 120001) law. Values above this line 
are due to numerical noise. 



/esc(^o) for reasonable widths of the emission line profile (typi- 
cally /esc varies by less than 0.2 dex for FWHM(LyQ') between 
50 and 200 km s '), since the resulting Hi absorption line pro- 
file is broader than FWHM for most column densities. In any 
case, we compute /sc assuming a Gaussian emission line profile 
with varying FWHM(LyQ') from 50 to 1000 km s"' for various 
applications. These values are provided in Table |2] 

4. 1 . Dependence of /esc on the physical parameters 

Since Lya photons can only be destroyed by absorption by dust 
particles, /esc mainly depends on the dust optical depth, de- 
scribed here by Tg. The predicted Lya escape fraction for all 
models is shown as a function of Eb-v in Fig. [T]for lines with 
input FWHM < 200 km s Here Eb-v is derived from the es- 
cape fraction in the continuum near Lya, assuming k,x - \2 and 
Rv - 4.05 from the ICalzetti et all (l2000h law. Very low escape 
fractions are obtained in some models. In such cases the line 
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Table 2. Derived quantities /esc, /es™'' ^^'^ Lye equivalent width for a constant input spectrum (flat continuum) (cols. 6-8) computed 
for all 6240 models (described by the parameters given in cols. 1-4) and for different FWHM values of the input Lya line (col 5). 
Note that the derived quantities are identical for cases of spherical infall with the same absolute radial velocity (|Vexpl). The full table 
is available only in the electronic version or on request from the author. 

h r„ logA'Hi Vexp FWHM(LyQ') U EW 

[kms-'] [cm--] [kms"'] [kms"'] [A] 



10. 


0.2 


16.0 


0. 


50. 


0.41473 


0.84103 


-0.13742 


10. 


0.2 


16.0 


0. 


100. 


0.57985 


0.84103 


-0.13742 


10. 


0.2 


16.0 


0. 


150. 


0.65799 


0.84103 


-0.13742 


10. 


0.2 


16.0 


0. 


200. 


0.70028 


0.84103 


-0.13742 
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Fig. 2. Predicted Lya escape fractions from our 6240 models as 
a function of A^hi- Colour codes for the symbols as in Fig.[T] 



profiles are dominated by a broad absorption (Voigt-like profile) 
and no distinctive trace of the intrinsic Lya emission line is de- 
tectable. For this reason we limit our plots, somewhat arbitrarily, 
at /esc ~ 10"^. This limit also corresponds to our numerical limit 
on the escape fraction per frequency bin. 

The Lya escape fraction has an upper limit of /esc ^ 
exp(-Ta) » exp(10 X Eb-v), corresponding to the pure atten- 
uation of the continuum. Below this value we see that /esc varies 
by several orders of magnitudes, when variations of the remain- 
ing parameters (Vgxp, A^hi, b) over a large range are allowed. In 
this case, multiple scattering effects of Lya on HI and on dust in- 
crease the probability of subsequent dust absorption of Lya pho- 
tons, reducing therefore /esc. For example, for a given Ta(EB-v), 
increasing the HI column density, A^hi, leads to wide range of 
/esc, allowing thus in particular low escape fractions. Similarly, 
a nearly static ISM (low expansion velocities Vexp) leads to lower 
Lya escape fractions, when all other parameters are the same. In 
both cases the increase of A^hi and the decrease of Vexp increases 
the number of Lya (and dust) scatterings due to radiation transfer 
effects, and hence the probability of absorption by dust for same 
radial dust optical depth. In particular we note that relatively low 



Lya escape fractions (/-sc ^ 0.1 or less) can even be obtained in 
situations with very little dust {Eb-v~ 0.02), provided the ISM 
has low velocities with respect to the UV source. 

In reality not all combinations of parameters may be reaUsed. 
For example, the dust-to-gas ratio 



[dust/ gas] = l0g(£'B_v/A^Hl) - log(£'B-v/A^Hl)Gal 

^ log(0.lTj-log(A^Hi) + 21.76 



(2) 



may be restricted within some range fro m the Galactic valu e 
log(A^Hi/£B-v)Gai = log(5.8 X lO^') cm^^ (|BohlinetalJll978h . 
As clear from Table [1] some of our models have extreme dust- 
to-gas ratios both above and below unity. In Figs. 1-4 we there- 
fore distinguish with different colors the models which dust- 
to-gas ratios within a factor 10 (2) of the Galactic value, i.e. 
with [dust/gas] 6 [-1, 1] ([-0.3,0.3]), and the remaining models. 
Among the models with dust-to-gas ratios close to the Galactic 
value (within a factor 2; i.e. among the red points) we have fur- 
ther highlighted in blue models with a typical outflow velocity 
of Vexp=200 km s~' to illustrate how e.g. velocity affects the Lya 
escape fraction. The remaining spread in /esc still shows how the 
other parameters (b, Nhi) affect the radiation transfer and hence 
the predicted Lya escape. 

4.2. Lya absorption 

Radiation transfer effects not only regulate the transmission of 
Lye line photons discussed above; they also affect the UV con- 
tinuum photons, "carving" thus broad Lya absorption lines of- 
ten observed in spectra of distant objects. Whereas scattering of 
photons out of the line-of-sight is principally at the origin of ab- 
sorption lines of the Lyor forest, true absorption by dust is the 
only processes capable of creating Lya absorption lines for the 
"closed" geometry and for integrated spectra we consider here 
(cf. Verhamme et al. 2006). 

In Fig. [3] we show the predicted equivalent width of Lya ab- 
sorption (EW < 0) from our model grid as a function of the 
Hi column density. Here EWiya has been computed assuming a 
constant (flat) UV spectrum and measuring EW over the interval 
[-6000,6000] km s ' covered by our modelfl Dust-free mod- 
els show EW - CB whereas otherwise \EW\ increases broadly 
with A^Hi, as expected. For a given column density, the strength 
of the Lya absorption increases with increasing dust-to-gas ratio 



- For the highest column densities, the Lya line is broader than this. 
Measuring EWiy„ over [-10000,10000] km s"' instead, increases EW 
by ~ 10-20 %. 

^ Few models show EWLya> but close to zero. This is due to nu- 
merical noise. 



Daniel Schaerer, Matthew Hayes, Anne Verhamme, Romain Teyssier: Lyo- radiation transfer models 



5 




Fig. 3. Predicted EWi,ya in absorption from our model grid as a 
function of A^hi (black points) for a flat (constant) input spec- 
trum, simulating a pure UV continuum. Colour codes for the 
symbols as in Fig. [1] Predictions at EWi,ya= correspond to 
dust free models. The dotted lines show EW^ya obtained from 



pure Voigt profiles for b 
right to left) 



10, 20, 40, 80, and 160 km (from 



(coded in color). It also depends on the other parameters aff^ect- 
ing the radiation transfer, here mostly Vgxp. 

For comparison, EW computed from Voigt profiles as func- 
tion of A^Hi and b (in the same fashion, i.e. over [-6000,6000] 
km s~') are shown by dotted lines in Fig. |3] As discussed by 
I Verhamme et al.l ( 2006h and shown e.g. from d etailed line pro- 
file fitting bv'Dessauges-Za vadskv et alJ (|2010'), radiation trans- 
fer models for the geometry adopted here predict in some cases 
weaker Lya absorption than expected for A^hi deduced by simple 
Voigt profiles. In other words, if applicable, our geometry could 
imply that A^hi measured from pure Voigt profile fits are underes- 
timated. However, the precise amount of this difference depends 
on the b value adopted for the Voigt fits, on the dust-to-gas-ratio, 
and on other parameters. 

Overall we note that the strength of the Lya absorption pre- 
dicted by our model grid covers well the range of observed 
EWiya values in Lyman Break Galaxies (LBGs), with equiva- 
lent widths down to ~ -20 to -30 A (e.g. lShaplev et al.ll2003h . 

As expected, there is no strong correlation between the es- 
cape fraction of Lya line photons and EWabs (see Fig. |4ll. The 
former measures the transmission of the photons close to line 
center, the latter is an overall measure of the net transmission 
of the UV flux. Therefore it is clear that cases with a strong ab- 
sorption line (EW^ya'S: 0) also show a low Lya escape fraction. 
However, for some models with low Hi column densities and 
high dust-to-gas ratios /esc can be quite low, but their overall ab- 
sorption remains weak (EW ~ -5 to A). These models, e.g., 
occupy the lower right corner of Fig. 2] Overall, models with a 
large EW in absorption correspond to high column densities (cf. 
Fig. [3]). As the Lya escape fraction mostly reflects the transmis- 
sion of photons close to line center, /esc decreases more rapidly 



Fig. 4. Predicted Lya escape fractions from our 6240 models as 
a function of the absorption equivalent width expected for a flat 
(constant) UV spectrum. Colour codes for the symbols as in Fig. 

than the equivalent width of the absorption line increases (in ab- 
solute terms). This explains the near vertical behaviour of the 
sequences of constant A^hi shown in Fig.|4] 

Since in reality, the intrinsic UV spectrum of star-forming 
galaxies is in gene ral composed of a UV contin uum plus Lya 
emission (see e.g. ISchaerer & Verhammel l2008i for synthetic 
spectra in this region). Fig. |4] schematically shows how the 
two components (line and continuum) are affected by radiation 
transfer effects. Approximately, the resulting emergent line is 
then the superposition of the absorption carved from the con- 
tinuum plus the remainder of the transmitted Lya emission 
line (with a positive equivalent width). To predict the detailed 
sha pe Lya profile, sh o wing a complex diversity as e.g . show n 



by Verhamme et al.1 (l2006l): ISc haerer & Ve rharnmel (l2008l) 



Dessauges-Zavads kv et all i2010 ): Vanzella et alj (l2010l) . flie 
full results from our radiation transfer simulations need to be 
used. We now briefly illustrate some profiles and variations with 
the model parameters. 

5. A library of theoretical Lya line profiles 

As already mentioned, the results from our radiation transfer 
models can be used to compute a posteriori the predicted spec- 
trum around Lya for any arbitrary input spectrum. To illustrate 
some Lya line profiles predicted by our extensive model grid, 
we show few selected results computed for a continuum plus 
Lya emission, assuming FWHM= 100 km s"' and EW^ya- 60 
A (except mentioned otherwise), and plotted for a spectral reso- 
lution of 150 km s"' . 

In Fig.|5]we show simulations for a rapidly expanding shell 
(Vexp = 300 km s"', assuming b - 40 km s"') with vary- 
ing Hi column densities (left panel), and varying dust content 
( right panel). Bo th panels are similar to Figs. 16 and 17 of 
iVerhamme et al.l (|g006), but we here show the predictions for 
a wider parameter space. Qualitatively the behaviour is as dis- 
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Ny, effects for v = 300 km/s, no dust dust effects for v = 300 km/s, log(N^,) = S0.2 cm^^ 




-1000 



1000 
V [km/s] 



2000 




-1000 



1000 
V [km/s] 



2000 



Fig. 5. Predicted normalised Lya line profiles for an expanding shell with Vexp - 300 km s ' , = 40 km s ' , and with varying Hi 
column densities (left panel), and varying dust content (right panel). All line profiles are normalised around ± 5900 km s Left: 
Shown are models with logA^m = 19 .3 (blue dotted), 19.9 (red dotted), 20.2 (black), 20.8 (green dashed) cm"^ comparable to those 
in Fig. 16 of IVerhamme et all (l2006l) . plus models with lower column densities (log A^hi - 16., 18., 19. cm"^, solid blue), and with 
higher lower column densities (log Ahi = 2 1 . 1 , 2 1 .4, 2 1 .7 cm"^, solid red). Right: Shown are models with t„ = 0. (black, top), 0. 1 
(blue dotted), 0.5 (red dashed), 1 . (green dashed) comparable to those in Fig. 16 of iVerhamme et al., (i2006) . plus models with higher 
optical depths (to - 1.5, 2., 3. (solid blue), and 4. (solid black)). 



cussed and explained by IVerhamme et al.l (l2006l) . In particu- 
lar, for models with different column densities, it is seen that 
the position of the most prominent peak of the profile shifts to 
^ (1.5 - 2) X Vexp for log Ahi ^ 10^" cm^^ (left panel). For the 
highest column densities the profile becomes very flat and broad, 
extending to very high velocities (> 1500 km s"'). The right 
panel illustrates the effect of dust (absorption) optical depths 
up to Ta - 4. While to first order the profile shape remains 
similar (when renormalised) but somew hat "sharpened" for rel- 
ativel y low optical depths (tq ^ 1. cf. IVerhamme et"ani2006i 
l2008HLaursen et a l." 2009b). the overall shape of the profile can 
be more significantly altered at higher optical depths. Here, e.g., 
the main peak gets nearly destroyed for Ta -4-.. 

Fig.|6]shows the line profile from an outflow with Vgxp - 150 
km s"', a relatively high column density (log Ahi = 20.8 cm"^), 
and for varying dust content (Tq = to 4.). Note how the in- 
creasing dust content progressively reduces the line to "carve" 
out a damped (Voigt-like) profile for the highest dust contents. 
Such damped profiles with some remaining Lyes' emission in the 
wings have e.g. been observed for relatively dusty LBGs such as 
cB58 and the 8 o'clock arc, which have successfully been mod- 
eled \yith simulations from this grid (c f. ISchaerer & Verhammg 
l2008tlDessauges-Zavadskv et alfcoiol) . 

As a final illustration of the large parameter space covered by 
our models, we show in Fig. |7] how the Lya profile varies with 
shell velocity in a very dusty, large column density environment. 
As expected, an emission line emerges out of the damped pro- 
file seen for a static shell (black) when the shell has a net radial 
flow velocity, and quite rapidly the profile becomes resembles a 
P-Cygni profile with blue shifted absorption and redshifted emis- 



sion, with the complex detailed shape whose origin has been dis- 
cussed in Verhamme et al. (2006). 

Again, note that our results also apply for symmetry reasons 
to models with spherical infall instead of outflows. In this case 
the predicted line profiles are simply inverted in velocity space. 

To make efficient use of this large spectral library we 
(MH) have also developed an automati c fitting tool , whic h 
has already been app l ied in iDess auges-Z avadskv et al.l (l2010l) : 
IVanzeUa etal] ( 1201 01) : iLidman et all ( 1201 ll) . Other appHcations 
will be presented elsewhere. 



6. Conclusions 

We have improved o ur 3D Lya and UV co ntinuum radiation 
transfer code McLya (Verha mme et al.ll2006l) and parallelised it 
for efficient use. Using this code we have computed a large grid 
of radiation transfer models for homogeneous spherically ex- 
panding (and by symmetry also infalling) shells containing Hi 
and dust surrounding a central source. With 6240 simulations, 
the grid covers a large parameter space with radial velocities 
from to 700 km s"', neutral hydrogen column densities A^hi 
from 10'^ to 10^''^ cm"^, Doppler parameters b - 10 to 160 
km s"', and dust optical depths Tg = 0. to 4 (corresponding ap- 
proximately to Eb-v~ 0.4). 

From the simulations we determine the predicted Lya escape 
fraction, the UV continuum attenuation, and we illustrate some 
of the Lya line profiles which are predicted by these models. 
The model predictions provide a useful basis for the interpreta- 
tion of Lya observations (line fluxes, escape fractions, equiva- 
lent widths, detailed line profiles and others) from star-forming 
galaxies, including Lyman alpha emitters (LAEs), Lyman break 
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Note how the increasing dust content progressively reduces the 
line to "carve" out a damped (Voigt-like) profile for the highest 
dust contents. 



galaxies (LBGs) and others. The models have already been suc- 
cess fuUy applie d and confronted to a variet y of observations (see 
e.g. lAteketal."2009a'; 'Hayes et al."2010',^20Tl'; 'Vanzellaetafl 
12010; Dessauges-Zavadskv et al., ,2010; .Garel et al.. 201 lb and 
are now made available publicly. 

Our models should also serve as an important base for com- 
parison for future models including other refinements (e.g. other 
geometries and velocity fields, non-homogeneous media) and 
further improvements. 
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